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Abstract — We study the problem of estimating vector- valued 
variables from noisy "relative" measurements. This problem 
arises in several sensor network applications. The measurement 
model can be expressed in terms of a graph, whose nodes 
correspond to the variables and edges to noisy measurements 
of the difference between two variables. We take an arbitrary 
variable as the reference and consider the optimal (minimum 
variance) linear unbiased estimate of the remaining variables. 

We investigate how the error in the optimal linear unbiased 
estimate of a node variable grows with the distance of the node 
to the reference node. We establish a classification of graphs, 
namely, dense or sparse in R'*, 1 < d < 3, that determines how 
the linear unbiased optimal estimation error of a node grows with 
its distance from the reference node. In particular, if a graph is 
dense in 1,2, or 3D, then a node variable's estimation error is 
upper bounded by a linear, logarithmic, or bounded function of 
distance from the reference, respectively. Corresponding lower 
bounds are obtained if the graph is sparse in 1, 2 and 3D. Our 
results also show that naive measures of graph density, such as 
node degree, are inadequate predictors of the estimation error. 
Being true for the optimal linear unbiased estimate, these scaling 
laws determine algorithm-independent limits on the estimation 
accuracy achievable in large graphs. 



I. Introduction 

Several applications in sensor and actuator networks lead 
to estimation problems where a number of variables are to be 
estimated from noisy measurements of the difference between 
certain pairs of them. Consider the problem of localization, 
where a sensor does not know its position in a global coor- 
dinate system, but can measure its position relative to a set 
of nearby nodes. These measurements can be obtained, for 
example, from range and angle data but are typically subjected 
to large noise (see Figure [TJ. In particular, two nearby sensors 
u and V located in a plane at positions Xu and Xy, respectively, 
have access to the measurement 
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where eu.v denotes measurement error The problem of interest 
is to use the C,u.v s to estimate the positions of all the nodes in 
a common coordinate system whose origin is fixed arbitrarily 
at one of the nodes. 

Similar estimation problems arise in time synchroniza- 
tion Ul 13, S and motion consensus in sensor-actuator net- 
works U; see Ed for an overview of these applications. 
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Fig. 1 . Relative position measurement in a Cartesian reference frame using 
range and angle measurements. A local compass at each sensor is needed to 
measure bearing with respect to a common North. Noisy measurements of 
the range ru,v and angle du,v between a pair of sensors u and v, which are 
denoted by fu,v and 9u,v, are converted to noisy measurements of relative 

position in the x — y plane as C,u,v = —lfu,v cos du,v , ru,v sin9u,v]'^ , 

with c = E[cos((59)], where Sg = 6 — 6 is the random error in the angle 
measurement. The division by c is needed to ensure that the noise in the 
measurement (^u,v is zero mean (see Appendix |ll}. The same procedure is 
performed for every pair of sensors that can measure their relative range and 
angle. The task then is to estimate the positions of all the nodes with respect 
to an arbitrary node in the network from the relative position measurements. 



Motivated by these applications, we study the problem of 
estimating vector valued variables from noisy measurements 
of the difference between them. In particular, denoting the 
variables of interest by {xi : i e V} where V := {1, 2, . . . }, 
we consider problems for which noisy relative measurements 
of the form ([TJ are available. The ordered pairs of indices 
{u, v) for which we have relative measurements form a set 
E that is a (typically strict) subset of the set 1^ x 1/" of all 
pairs of indices. Just with relative measurements, the Xu's can 
be determined only up to an additive constant. To avoid this 
ambiguity, we assume that a particular variable (say Xq) is 
used as the reference, which is therefore assumed known. The 
problem of interest is to estimate the remaining node variables 
from all the available measurements. 

The measurement equations ([T]i can be naturally associated 
with a directed graph Q = {'V ^T, ) with an edge from node u 
to V if the measurement (^u.v is available. The graph Q is called 
the measurement graph, and each vector Xu, u G is called 
the u-th node variable. The measurement noise ee,e G £ is 
assumed zero mean and spatially uncorrelated, i.e., E[ee] — 
Ve e £ and ^[ee^] = if e ^ e. 

In this paper we investigate how the structure of the graph Q 
affects the quality of the optimal linear unbiased estimate Xu 
of Xu, measured in terms of the covariance of the estimation 
error E„^o := E[(a;K — Xu){xu — Xu)^]- The optimal linear 
unbiased estimate refers to the one obtained with the classical 
best linear unbiased estimator (BLUE), which achieves the 
minimum variance among all linear unbiased estimators i^. 



We examine the growth of the BLUE error variance of a node 
u as a function of its distance to the reference node. 

We are interested in the growth of error with distance in 
large graphs, for which infinite graphs (with a countably infi- 
nite number of nodes and edges) serve as proxies. This paper 
is focused on infinite graphs because the absence of boundary 
conditions in infinite graphs allows for more complete and 
simpler results. Using infinite graphs as proxies for large finite 
graphs is theoretically justified by the fact that the BLUE 
error variance of a node variable Xu in a large but finite 
subgraph of an infinite graph is arbitrarily close to the BLUE 
estimation error in the infinite graph, as long as the finite graph 
is sufficiently large. This convergence result was established 
in 0]. 

When the measurement graph is a tree, there is a single path 
between the u* node and the reference node and one can show 
that the covariance matrix of the estimation error is the sum 
of the covariance matrices associated with this path. Thus, for 
trees, the variance of the BLUE estimation error of Xu grows 
linearly with the distance from node u to the reference node. 
It turns out that for graphs "denser" than trees, with multiple 
paths between pairs of nodes, the variance of the optimal linear 
unbiased estimation error can grow slower than hnearly with 
distance. 

In this paper, we introduce a novel notion of denseness for 
graphs that is needed to characterize how the estimation error 
grows with distance. In classical graph-theoretic terminology, 
a graph with n vertices is called dense if its average node 
degree is of order n, and is called sparse if its average node 
degree is a constant independent of n |8]. We recall that the 
degree of a node is the number of edges incident on it (an 
edge (u, v) is said to be incident on the nodes u and v). Other 
notions of denseness include geo-denseness introduced by 
which requires uniform node density (nodes per unit area) but 
does not consider the role of edges. Accuracy of localization 
from distance-only measurements have been extensively stud- 
ied in the sensor networks literature, ty pica lly by evaluating 
the Cramer-Rao lower bound (see iH llil El El Q and 
references therein). In many of these studies, graph density 
(as measured by node degree or node density) is recognized to 
affect estimation accuracy However, we will see 

through examples in Remark |2] that for the estimation problem 
considered in this paper, such notions of denseness are not 
sufficient to characterize how the estimation error grows with 
distance. 

A key contribution of this paper is the development of 
suitable notions of graph denseness and sparseness that are 
useful in determining BLUE error scaling laws. These notions 
exploit the relationship between the measurement graph and 
a lattice. We recall that the rf-dimensional square lattice 
is defined as a graph with a node in every point in with 
integer coordinates and an edge between every pair of nodes 
at an Euclidean distance of 1 (see Figure |4] for examples). 
The error scaling laws for a lattice measurement graph can be 
determined analytically by exploiting symmetry. It turns out 
that when the graph is not a lattice, it can still be compared 



to a lattice. Intuitively, if after some bounded perturbation in 
its node and edge set, the graph looks approximately like a 
rf-dimensional lattice, then the graph inherits the denseness 
properties of the lattice. In that case the error covariance for 
the lattice can still be used to bound the error covariance in 
the original graph. 

Our classification of dense and sparse graphs in W^, d E 
{1,2,3}, characterizes BLUE error scaUng laws. For dense 
graphs, they provide upper bounds on the growth rate of the 
error, while for sparse graphs, they provide lower bounds. The 
precise growth rates depend on which dimension the graph is 
dense or sparse in. When a graph is dense in ID, 2D, or 3D, 
respectively, the error covariance of a node is upper bounded 
by a linear, logarithmic, or bounded function, respectively, of 
its distance from the reference. On the other hand, when a 
graph is sparse in ID, 2D, or 3D, the error covariance of a 
node is lower bounded by a linear, logarithmic, and bounded 
function, respectively, of its distance from the reference. Our 
sparse graphs are also known as "graphs that can be drawn in 
a civilized manner" according to the terminology introduced 
by Doyle and Snell lll5ll in connection with random walks. 

The BLUE error scaling laws derived in this paper provide 
an algorithm-independent limit to the estimation accuracy 
achievable in large networks, since no linear unbiased esti- 
mation algorithm can achieve higher accuracy than the BLUE 
estimator For example, when a graph is sparse in ID, the 
BLUE estimation error covariance grows at least linearly 
with the distance from the reference. Therefore the estimation 
accuracy will be necessarily poor in large ID sparse graphs. 
On the other hand, when a graph is dense in 3D, the BLUE 
estimation error of every node variable remains below a 
constant, even for nodes that are arbitrarily far away from 
the reference. So accurate estimation is possible in very large 
3D dense graphs. 

The results in this paper are useful for the design and 
deployment of ad-hoc and sensor networks. Since we now 
know what structural properties are beneficial for accurate 
estimation, we can strive to achieve those structures when 
deploying a network. Specifically, we should try to achieve 
a dense-in-R'* structure, with d as large as possible, for high 
accuracy estimation. Since the scaling laws are true for the 
optimal linear unbiased estimator, they can also help designers 
determine if design requirements are achievable. For example, 
if the requirement is that the estimation accuracy should not 
decrease with size, no matter how large a network is, the 
network must be dense in R'', d > 3 for such a requirement 
to be satisfied. 

Our results also expose certain misconceptions that exist in 
the sensor network literature about the relationship between 
graph structure and estimation error In Section III-BI we 
provide examples that expose the inadequacy of the usual mea- 
sures of graph denseness, such as node degree, in determining 
scaling laws of the estimation error. 

In practice, more than one reference node (commonly 
referred to as anchors) may be used. We only consider the 
case of a single reference node since scaling laws with a single 
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reference provide information on how many reference nodes 
may be needed. For example, since the estimation error in a 
3D dense graph is bounded by a constant, one reference node 
may be enough for such a graph. 

While we do not discuss the computation of the optimal 
linear unbiased estimates in this paper, we have developed dis- 
tributed algorithms to compute these estimates with arbitrary 
precision (see (|5|] and references therein). These algorithms 
are distributed in the sense that every node computes its 
own estimate and the information needed to carry out this 
computation is obtained by communication with its neighbors. 

A preliminary version of some of the results in this paper 
was presented in 116}. However, iflil] used stricter assumptions 
to establish the upper bounds on error growth rates. Moreover, 
only sufficient conditions were obtained in 1 16] for some of the 
error scaling laws to hold; whereas here we derive necessary 
and sufficient conditions. 

Organization: The rest of the paper is organized as follows. 
Section HI] describes the problem and summarizes the main 
results of the paper Section describes key properties of 
dense and sparse graphs. Section |IV] briefly describes the 
analogy between BLUE and generalized electrical networks 
from |7] that is needed to prove the main results. Section IVl 
contains the proof of the main result of the paper. Section |VI] 
deals with the question of how to check if a graph possesses 
the denseness/sparseness properties. The paper ends with a 
a few final conclusions and directions for future research in 
Section [Vn] 

II. Problem Statement and Main Results 

Recall that we are interested in estimating vector-valued 
variables G M*^, u E 'V := {1, 2, . . . }, from noisy relative 
measurements of the form: 

where eu,v denotes a zero-mean measurement noise and "E is 
the set of ordered pairs {u,v) for which relative measurements 
are available. The node set 1^ is either finite, or infinite but 
countable. We assume that the value of a particular reference 
variable Xq is known and without loss of generality we take 
Xo = 0. The node set 'V and the edge set £ together define 
a directed measurement graph Q = , £ ). 

The accuracy of a node variable's estimate, measured in 
terms of the covariance of the estimation error, depends on 
the graph Q as well as the measurement errors. The covariance 
matrix of the error „ in the measurement Qu.v is denoted by 
Pu,v, i-e., Pu,v '■= E[eu,Deu i)]- We assume that the measure- 
ment errors on different edges are uncorrected, i.e., for every 
pair of distinct edges e, e € £, E[eeeJ] = 0. The estimation 
problem is now formulated in terms of a network {Q, P) 
where P : ^ is a function that assigns to each edge 
{u, w) e £ the covariance matrix P^^v of the measurement 
error associated with the edge {u, v) in the measurement graph 
Q. The symbol S'^^ denotes the set of kxk symmetric positive- 
definite matrices. 



As discussed in Section U our results are stated for infinite 
networks. The following conditions are needed to make sure 
that the estimation problem is well posed and that the estimates 
satisfy appropriate convergence properties to be discussed 
shortly: 

Assumption 1 (measurement network): The measurement 
network (Q, P) satisfies the following properties: 

1) The graph Q is weakly connected, i.e., it is possible to go 
from every node to every other node traversing the graph 
edges without regard to edge direction. 

2) The graph Q has a finite maximum node degre^H. 

3) The edge-covariance function P is uniformly bounded, 
i.e., there exists constant symmetric positive matrices 

-Pmin, -Pmax SUCh tiiat Pmin < Pe < ^'max, Vc G £ . □ 

In the above, for two matrices A, B € R''^'', A > B (A > 
B) means A — i? is positive definite (semidefinite). We write 

A < B {A < B) if -A> -B {-A > -B). 

We also assume throughout the paper that measurement 
graphs do not have parallel edges. A number of edges are 
said to be parallel if all of them are incident on the same pair 
of nodes. The condition of not having parallel edges is not 
restrictive since parallel measurements can be combined into 
a single measurement with an appropriate covariance, while 
preserving the BLUE error covariances (see Remark [3]). 

Given a finite measurement network [G&nitcT P), where 
^finite contains the nodes u and o, it is straightforward to 
compute the BLUE estimate i„(5finito) of the unknown node 
variable a;„ in the network {Qfinite, P), as described in and 
the covariance matrix of the estimation error SM,o(^finitc) := 
E[(a;„ — XyMxu — Xu)'^] exists as long as t/finito is weakly 
connected ||7|]. Due to the optimality of the BLU estimator, 
SM,o(^finite) is the minimum possible estimation error covari- 
ance that is achievable by any linear unbiased estimator using 
all the measurements in the graph ^/finite- 

When the measurement graph is infinite, the BLUE error 
covariance S]„.o for a node variable Xu is defined as 

Su,o = inf S„.o(tJfinito), (3) 

Sfinitc 

where the infimum is taken over all finite subgraphs ^/finite of 
Q that contain the nodes u and o. We define a matrix M to 
be the infimum of the matrix set S C S*^^, and denote it by 

M = inf 5, (4) 

if M < A for every matrix A E S, and for every positive 
real e, there exists a matrix B E S such that M + elk > B. 
Under Assumption [T] it was shown in |7] that the infimum 
in ^ always exists. In this case, ^ means that the BLUE 
covariance T^u.o is the the lowest error covariance that can be 
achieved by using all the available measurements. 

In the sequel, we determine how the BLUE covariance I]„ „ 
grows as a function of the distance of node u to the reference 
o, and how this scaling law depends on the structure of the 

' The degree of a node is the number of edges that are incident on the node. 
An edge {u,v) is said to be incident on the nodes u and v. 
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measurement graph G- To this effect we start by providing a 
classification of graphs that is needed to characterize the error 
scahng laws. 

A. Graph Denseness and Sparseness 

We start by introducing graph drawings, which will later 
allow us to define dense and sparse graphs. 

1 ) Graph Drawings: The drawing of a graph Q = ( 1^ , £ ) 
in a d-dimensional Euclidean space is obtained by mapping 
the nodes into points in M.'^ by a drawing function f : "V ^ 
W'-. A drawing is also called a representation [17] or an 
embedding [8]. For a particular drawing /, given two nodes 
u,v ^ 'P the Euclidean distance between u and v induced by 
the drawing / : 1/" ^ M"^ is defined by 

df{u,v) ■.= \\f{v)^f{u)\\, 

where || • || denoted the usual Euclidean norm in M.'^. It is 
important to emphasize that the definition of drawing allows 
edges to intersect and therefore every graph has a drawing 
in every Euclidean space. In fact, every graph has an infinite 
number of drawings in every Euclidean space. However, a 
particular drawing is useful only if it clarifies the relationship 
between the graph and the Euclidean space in which it is 
drawn. In what follows, given two nodes u and v, dg{u,v) 
denotes the graphical distance between u and v, i.e., the 
number of edges in the shortest path between u and v. The 
graphical distance dg is evaluated without regards to edge 
directions, which are immaterial in determining BLUE error 
covariances (see Remark H]). 

For a particular drawing / and induced Euclidean distance 
d/ of a graph Q — four parameters are needed to 

characterize graph denseness and sparseness. The minimum 
node distance, denoted by s, is defined as the minimum 
Euclidean distance between the drawing of two nodes 

s inf df{u, v). 

The maximum connected range, denoted by r, is defined as 
the Euclidean length of the drawing of the longest edge 

r := sup df{u, v). 
(u,i>)eE 

The maximum uncovered diameter, denoted by 7, is defined 
as the diameter of the largest open ball that can be placed in 
M.'^ such that it does not enclose the drawing of any node 



s = 1 
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sup|(5 : 3Bs s.t. f{u) ^ Bs^u e -J^j, 



where the existential quantification spans over the balls Bs in 
M.'^ with diameter S and centered at arbitrary points. Finally, 
the asymptotic distance ratio, denoted by p, is defined as 

p := lim inf / ^/(^' ^ 1/ and dg{u,v) > nX. 

n^oo ldg[U,V) J 

Essentially p provides a lower bound for the ratio between 
the Euclidean and the graphical distance for nodes that are far 
apart. The asymptotic distance ratio can be thought of as an 




Fig. 2. A drawing of a graph in 2D Euclidean space, and the corresponding 
denseness and sparseness parameters. Since the minimal distance between 
any two nodes is 1, the minimum node distance is s = 1. Since the longest 
edge is between u* and v* , the maximum connected range is r = vTO. The 
diameter of the largest ball that can fit inside the drawing without enclosing 
any node is 2. so the maximum uncovered diameter is thus 7 = 2. The 
minimal ratio between the Euclidean and graphical distance of a pair of nodes 
is achieved by the pair p* ,q* , hence the asymptotic distance ratio is p = 
df(p'',q*)/dgip*,q*) = 1/5. 



inverse of the stretch for geometric graphs, which is a well- 
studied concept for finite graphs OlSll . 



The two parameters p and r defined above are especially 
useful to compare graphical and Euclidean distances, as stated 
in the following result. 

Lemma 1 (Euclidean vs. graphical distances): The follow- 
ing two statements are equivalent: 

1) The asymptotic distance ratio p is strictly positive. 

2) There exist constants a > 0, /? > for which 



dg{u,v) < adf{u,v) + (3, \fu,vG'P'. 

Similarly, the following statements are equivalent: 

1) The maximum connected range r is finite. 

2) There exist constants a > 0, f3 > for which 
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□ 



df{u,v)<adg{u,v)+P, Wu,v£1/. 
The proof of this lemma is provided in Appendix U 
2) Dense and Sparse Graphs: We call the drawing of a 
graph with finite maximum uncovered diameter (7 < 00) and 
positive asymptotic distance ratio (p > 0) a dense drawing. We 
say that a graph Q is dense in M'' if there exists a dense drawing 
of the graph in R^. Graph drawings for which the minimum 
node distance is positive (s > 0) and the maximum connected 
range is finite (r < 00) are called civilized drawings lllSil . A 
graph Q is said to be sparse in M'' if there exists a civilized 
drawing in M''. 

It follows from these definitions and Lemma \T\ that if a 
graph is dense in R'', then it has enough nodes and edges 
so that it is possible to draw it in R'' in such a way that 
its nodes cover R'' without leaving large holes (finite 7), 
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and yet a small Euclidean distance between two nodes in the 
drawing guarantees a small graphical distance between them 
(positive p, which implies (|5]l). On the other hand, if a graph 
that is sparse in W^, then one can draw it in so as to 
keep a certain minimum separation between nodes (positive 
s) without making the edges arbitrarily long (finite r). It also 
follows from the definitions that a graph must be infinite to be 
dense in any dimension, and a finite graph is sparse in every 
dimension. 

A graph can be both dense and sparse in the same dimen- 
sion. For example, the d-dimensional lattice is both sparse and 
dense in M''. However, there is no civilized drawing of the d- 
dimensional lattice in for any d' < d. Moreover, there 
is no dense drawing of the d-dimensional lattice in M.'^ for 
every d > d. This means, for example, that the 3D lattice in 
not sparse in 2D and is not dense in 4D. In general, a graph 
being dense in a particular dimension puts a restriction on 
which dimensions it can be sparse in. The next result, proved 
in Section |Vl] states this precisely. 

Lemma 2: A graph that is dense in R'' for some d > 2, 
cannot be sparse in R'' for every d' < d. □ 

Remark 1 (historical note): In the terminology of Doyle 
and Snell |15t]. sparse graphs (as defined here) are said to 
be graphs "that can be drawn in a civilized manner". In this 
paper we refer to such graphs as sparse graphs since they are 
the antitheses of dense graphs. □ 



(a) A 3-fuzz of a ID lattice 




(b) A triangular lattice 




(c) A 3D lattice 



Fig. 3. Three measurement graphs that show vastly different scaling laws of 
the estimation eiTor, whereas each has the same node degree for every node. 
Furthermore, they are all "sparse" according to traditional graph-theoretic 
terminology (see the discussion on graph denseness in Section 



B. Error Scaling Laws 

The concepts of dense and sparse graphs allow one to 
characterize precisely how the BLUE error covariance o 
grows with the distance from the node u to the reference o. 
The next theorem, which establishes the BLUE error scaling 
laws for dense and sparse graphs, is the main result of the 
paper The proof of the theorem is provided in Section FVl 

Before we present the theorem, we need to introduce some 
notation. The asymptotic notations ri( ) and O(-) are used for 
matrix valued functions in the following way. For a matrix- 
valued function 9 : M ^ R'^ ^ and a scalar- valued function p : 
R ^ R, the notation g{x) — 0{p{x)) means that there exists 
a positive constant Xo and a constant matrix A G S''^ such 
that g{x) < Ap{x) for all x > Xo- Similarly, g{x) = il{p{x)) 
means there exists a positive constant Xq and a constant matrix 
B e such that g{x) > Bp{x) for all x > Xo- Recall that 
S''^ is the set of all fc x A; symmetric positive definite matrices. 

Theorem 1 (Error Scaling Laws): Consider a measurement 
graph Q — that satisfies Assumption [T] with a 

reference node o e 1^. The BLUE error covariance o for 
a node u obeys the scaling laws shown in Table U □ 

A graph can be both sparse and dense in a particular 
dimension, in which case the asymptotic upper and lower 
bounds are the same. For a graph that is both sparse and dense 
in R"*, the error covariance grows with distance in the same 
rate as it does in the corresponding lattice Z^. 

Remark 2 (Counterexamples to conventional wisdom): As 
noted in Section U the average node degree of a graph or 



the number of nodes and edges per unit area of a deployed 
network are often used as measures of graph denseness. 
However, these measures do not predict error scaling laws. 
The three graphs in Figure |3] offer an example of the 
inadequacy of node degree as a measure of denseness. This 
figure shows a 3-fuzz of the ID lattice (see Section |lll] for 
the formal definition of a fuzz), a triangular lattice, and a 
3-dimensional lattice. It can be verified from the definitions 
in Section III-A.2I that the 3-fuzz of the ID lattice is both 
dense and sparse in R, the triangular lattice is dense and 
sparse in R^, and the 3D lattice is dense and sparse in R'^. 
Thus, it follows from Theorem [T] that the BLU estimation 
error scales linearly with distance in the 3-fuzz of the ID 
lattice, logarithmically with distance in the triangular lattice, 
and is uniformly bounded with respect to distance in the 3D 
lattice, even though every node in each of these graphs has 
the same degree, namely six. □ 
We note that the notion of geo-denseness introduced in ^ 
is also not useful for characterizing error scaling laws since 
geo-denseness considers node density alone without regard to 
the edges. 

III. Dense and Sparse Graphs 

This section establishes an embedding relationship between 
dense and sparse graphs and lattices, which is needed to prove 
Theorem [1] Roughly speaking, a graph Q can be embedded 
in another graph ^ if ^ contains all the nodes and edges of 
Q, and perhaps a few more. The usefulness of embedding in 
answering the error scaling question is that when Q can be 
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TABLE I 

COVARIANCE T,u,o OF Xu' S BLUE ESTIMATION ERROR FOR GRAPHS THAT ARE DENSE OR SPARSE IN R''. IN THE TABLE, df^{u, o) DENOTES THE 

Euclidean distance between node u and the reference node o induced by a drawing : V — ♦ R'* that establishes the graph's 

DENSENESS IN THE EUCLIDEAN SPACER'', AND dji^{u, o) DENOTES THE EUCLIDEAN DISTANCE INDUCED BY A DRAWING : 1/" ^ R'' THAT 

ESTABLISHES THE GRAPH'S SPARSENESS. 



Euclidean space 
and graph example 


Covariance matrix S^.o of the estimation error of 
Xu in a sparse graph with a sparse drawing 


Covariance matrix o of the estimation error of 
Xu in a dense graph with a dense drawing 






— R 

- R2 

1 R3 




S„,o(e) =o(logd/,(w,o)) 



embedded in Q, the BLUE error covariances in Q are larger 
than the corresponding ones in Q (this statement will be made 
precise in Theorem |5] of Section HVl i. 

The h-fuzz of a graph Q, introduced by Doyle and Snell 
1 1511 . is a graph with the same set of nodes as Q but with a 
larger set of edges. Specifically, given a graph Q and a positive 
integer h, a /i-fuzz of Q, denoted by Q'^'^\ is a graph that has 
an edge between two nodes u and v whenever the graphical 
distance between these nodes in Q is less than or equal to h. 

We say that a graph Q = can be embedded in 

another graph Q ~ {'P','E) if C I', and, whenever there 
is an edge between two nodes in Q, there is an edge between 
them in Q. More precisely, Q can be embedded in Q if there 
exists an injective map 77 : — > 1^ such that for every 
{u,v) £ £, either {ri{u) , ri{v)) e £ or {ri{v) , ri{u)) £ £. In 
the sequel, we use G C G denote that G can be embedded 
in G- 

A. Relationship with lattices and Euclidean spaces 

The next theorems (Theorem |2] and |3]l show that sparse 
graphs can be embedded in fuzzes of Lattices, and fuzzes of 
dense graphs can embed lattices. In these two theorems we 
use <izd(') to denote the graphical distance in the lattice 
and df{-) to denote the Euclidean distance in K*^ induced by 
the drawing /. 

Theorem 2 (Sparse Embedding): A graph G — {"^ is 
sparse in R'' if and only if there exists a positive integer h such 
that G C Z[,''\ Moreover, if / : 1^ -> R"^ is a civilized drawing 
of G in R'', then there exists an embedding 77 : 'V ^ I'za so 
that Vw, V e 1^, 



dzAv{u),v{v)) > Vd(-df{u,v) - 2 



(6) 



where s is the minimum node distance in the /-drawing of G- 

□ 

In words, the theorem states that G is sparse in if and 
only if G can be embedded in an h-fnzz of a d-dimensional 
lattice. The significance of the additional condition (|6]l is that 
if the Euclidean distance between a pair of nodes u and w in a 



civilized drawing of the graph is large, the graphical distance 
in the lattice between the nodes that correspond to u and v 
must also be large. 

The first statement of Theorem |2] is essentially taken 
from [fisll, where it was proved that if a graph can be drawn 
in a civilized manner in R'^, then it can be embedded in a 
/i-fuzz of a d-lattice, where h depends only on s and r. A 
careful examination of the proof in |il5J reveals that it is not 
only sufficient but also a necessary condition for embedding in 
lattice fuzzes. The proof of this theorem is therefore omitted. 

Theorem 3 (Dense Embedding): A graph G = is 
dense in R'' if and only if there exists finite, positive integers 
h and c such that the following conditions are satisfied 

(i) ^C'' D Zd, and, 

(ii) if 77 : T/zd ^ is an embedding of Z^ into G^^K then, 
Vm e 1^, 3u e vi'^Za) ^ 1^ such that dg{u,u) < c. 

Moreover, if / : ^ R'' is a dense drawing of G in R'', 
then the embedding function ij in (ii) can be chosen so that 

Vm, u e 1^, we can find Uz,Vz € 'V'za satisfying 



dg{u,r]{uz)) < c, dg{v,r]{vz)) < c 
dza{uz,Vz) <4:d+ ^df{u,v) 



(7) 



where 7 is the maximum uncovered diameter of the /-drawing 
of □ 
In words, the two conditions state that G is dense in R"^' if 
and only if (i) the d-dimensional lattice can be embedded in 
an ft,-fuzz of G for some positive integer h and (ii) every node 
of G that is not the image of a node in Z^ is at a uniformly 
bounded graphical distance from a node that is the image of 
a node in Zd- The significance of (|7| is that not only we can 
find for every node in a close-by node that has a pre-image 
in the lattice, but also these close-by nodes can be so chosen 
so that if the Euclidean distance between a pair of nodes u 
and V in the drawing is small, then the graphical distance in 
the lattice between the pre-images of their close-by nodes is 
small as well. 
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IV. Electrical Analogy 

A crucial step in proving the main results of this paper is 
the analogy introduced in |i7J between the BLU estimation 
problem and an abstract electrical network, where currents, 
potentials and resistances are matrix valued. 

A generalized electrical network (t/, R) consists of a graph 
Q = (V, E) (finite or infinite) together with a function R : 
£ that assigns to each edge e € £ a symmetric 

positive definite matrix R^ called the generalized resistance 
of the edge. 

A generalized flow from node u €z I' to node u G 1^ with 
intensity j G M'''^''' is an edge-function j : £ ^ M'''^''' such 
that 



u 



(p,9)6E 
P=P 



E 

(?,p)eE 



h,p 







P 

p^v Vpe 'U. (8) 

otherwise 



A flow j is said to have finite support if it is zero on all but a 
finite number of edges. We say that a flow i is a generalized 
current when there exists a node-function y : V ^ R*^^*^ for 
which 



R. 



'U,V''U,V 



V(u,w)G£. 



(9) 



The node-function V is called a generalized potential as- 
sociated with the current i. Eq. dS) should be viewed as a 
generalized version of Kirchhoff's current law and can be 
interpreted as: the net flow out of each node other than u 
and V is equal to zero, whereas the net flow out of u is equal 
to the net flow into v and both are equal to the flow intensity j. 
Eq. (|9]l provides in a combined manner, a generalized version 
of Kirchhoff's loop law, which states that the net potential drop 
along a circuit must be zero, and Ohm's law, which states that 
the potential drop across an edge must be equal to the product 
of its resistance and the current flowing through it. A circuit is 
an undirected path that starts and ends at the same node. For 
fc = 1, generalized electrical networks are the usual electrical 
networks with scalar currents, potentials, and resistors. 

A. Effective Resistance and BLUE Error Covariance 

It was shown in 101 that when a current of intensity i G 
jgfcxfc flows from node u to node v, the resulting generalized 
current j is a linear function of the intensity i and there exists 
a matrix 

i??*, e such that 



neff 



ViG 



c fe X fc 



(10) 



We call the matrix the generalized effective resistance 

between u and v. In view of this definition, the effective 
resistance between two nodes is the generalized potential 
difference between them when a current with intensity equal 
to the identity matrix is injected at one node and extracted 
at the other, which is analogous to the definition of effective 
resistance in scalar networks ifisll . Note that the effective 
resistance between two arbitrary nodes in a generalized net- 
work is a symmetric positive definite matrix as long as the 
network satisfies Assumption [T] whether the network is finite 
or infinite iQl- 



Generalized electrical networks are useful in studying the 
BLU estimation error in large networks because of the follow- 
ing analogy between the BLU estimation error covariance and 
the generalized effective resistance. 

Theorem 4 (Electrical Analogy, from [f]): Consider a 
measurement network {Q, P) satisfying Assumption [T] with 
Q = ( 1^ , E ) and a single reference node o G 1^ . Then, for 
every node u G 1^ \ {o}, the BLUE error covariance Y,u,o 
defined in Q is a symmetric positive definite matrix equal to 
the generalized effective resistance Rf/^ between u and o in 
the generaUzed electrical network {Q,P): 



ncff 
Ru,o- 



□ 



Remark 3: In an electrical network, parallel resistors can 
be combined into one resistor by using the parallel resistance 
formula so that the effective resistance between every pair of 
nodes in the network remain unchanged. The same can be done 
in generalized electrical networks |19]. The analogy between 
BLUE covariance and effective resistance means that parallel 
measurement edges with possibly distinct measurement error 
covariances can be replaced by a single edge with an equiva- 
lent error covariance, so that the BLUE error covariances of all 
nodes remain unchanged. This explains why the assumption 
of not having parallel edges made at the beginning is not 
restrictive in any way. 

B. Graph Embedding and Partial Ordering of BLUE Covari- 
ances 

Effective resistance in scalar electrical networks satisfies 
Rayleigh's Monotonicity Law, which states that the effective 
resistance between any two nodes can only increase^ if the 
resistance on any edge is increased, and vice versa rfisll . The 
next result (proved in [7]), states that the same is true for 
generalized networks, whether finite or infinite. 

Theorem 5 (Rayleigh's Monotonicity Law Consider 
two generalized electrical networks {G,R) and {G,R) with 
graphs Q = {']/,'£) and Q = ('J^,£), respectively, such that 
both the networks satisfy Assumption [T] Assume that 

1) Q can be embedded in Q, i.e., Q C G, and 

2) Re > Re for every edge e G £ . 

Then, for every pair of nodes u,v £ 1^ of Q, 



off 



where R'^^ and R'^^ are the effective resistance between u 
and V in the networks {G, R) and {Q, R), respectively. □ 

The usefulness of Rayleigh's Monotonicity Law in an- 
swering the error scaling question becomes apparent when 
combined with the Electrical Analogy. It shows that when Q 
can be embedded in Q, the BLUE error covariances in Q are 
lower bounded by the error covariances in Q. Intuitively, since 
Q has only a subset of the measurements in Q, the estimates 
in Q are less accurate than those in Q. 

Remark 4: Although the graph Q that defines the electrical 
network {Q, R) is directed, the edge directions are irrelevant 
in determining effective resistances. This is why Rayleigh's 
Monotonicity Law holds with graph embedding, which is 
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(a) ID lattice Zi (b) 2D lattice Z2 (c) 3D lattice Z3 

Fig. 4. Lattices. 



insensitive to edge directions. The electrical analogy also 
explains why the edge directions are irrelevant in determining 
error covariances. □ 

C. Triangle Inequality 

Matrix-valued effective resistances satisfy a triangle in- 
equality, which will be useful in proving the error scaling 
laws in Section |V] It is known that scalar effective resistance 
obeys triangle inequality, and is therefore also referred to as 
the "resistance distance" [20]. Although the result in ll20ll 
was proved only for finite networks, it is not hard to extend 
it to infinite networks. The following simple extension of 
the triangle inequality to generalized networks with constant 
resistances on every edge was derived in ifioll : 

Lemma 3 (Triangle Inequality): Let (Q, Ro) be a general- 
ized electrical network satisfying Assumption [T] with a con- 
stant resistance Ro E S'^^ on every edge of Q. Then, for every 
triple of nodes u,v,w in the network. 



R 



lOff 



□ 



D. Effective Resistances in Lattices and Fuzzes 

Recall that given a graph Q and a positive integer h, the 
/i-fuzz of Q, denoted by Q'-'^K is a graph that has an edge 
between two nodes u and v whenever the graphical distance 
between them in Q is less than or equal to h. 

An /i-fuzz will clearly have lower effective resistance than 
the original graph because of Rayleigh's Monotonicity Law, 
but it is lower only by a constant factor as stated in the 
following result, which is a straightforward extension to the 
generalized case of a result about scalar effective resistance 
established by Doyle and Snell (see the Theorem on page 103, 
as well as Exercise 2.4.9, in BlSll ). The interested reader can 
find a proof in ifioll . 

Lemma 4: Let {Q, Ro) be a generalized electrical network 
satisfying Assumption[T]with a constant generalized resistance 
Ro G on its every edge. Let {G'^^\Ro) be the electrical 
network similarly constructed on G'^^\ the /i-fuzz of G. For 
every pair of nodes u and v in 1/ , 

aRZiQ)<KUs^'^)<KliQ), 

where is the effective resistance in the network (•, -Ro) 

and a e (0, 1] is a positive constant that does not depend on 
u and V. □ 



The following lemma establishes effective resistances in d- 
dimensional lattices and their fuzzes. 

Lemma 5: For a given positive integer h, consider the 
electrical network {Z^J^\ Ro) with a constant generalized 
resistance Ro G S'^^ at every edge of the /i-fuzz of the 
d-dimensional square lattice Z^. The generalized effective 
resistance R'^y between two nodes u and v in the electrical 

network {t}'^\Ro) satisfies 

1) K«(z(''') = e(dz,K^;)) 

2) i?^«(z(''') = e(logdz,(u,«)), 

3) K«(z('')) = e(i). □ 

Proof of Lemma \5l The scalar effective resistance in ID , 
2D, and 3D lattices follow linear, logarithmic and bounded 
growth rates, respectively | 2ll 12211 . Using these results, it was 
established in [7] that the matrix effective resistances in these 
lattices have the same scaling laws (see Lemma 5 in Thus, 
ID, 2D, and 3D lattices with matrix-valued resistances have 
linear, logarithmic, and bounded scaling laws for the effective 
resistance, which is the result with h ~ 1. The case h > 1 
follows from the application of Lemma |4] 

The slowing down of the growth of the effective resistance 
as the dimension increases can be attributed to the fact that the 
number of paths between each pair of nodes is larger in higher 
dimensional lattices. The scaling laws for effective resistance 
in lattices and their fuzzes also have intimate connections to 
the change from recurrence to transience of random walks in 
lattices as the dimension changes from d = 1,2 to d > S jisll . 



V. Proof of Theorem[T] 

We now prove Theorem[T]by using the tools that have been 
developed so far. The following terminology is needed for the 
proofs. For functions g : R ^ R'^^'^ and p : M ^ M, the 
notation g{y) — Q{p{y)) means that g{y) — n{p{y)) and 
g{y) — 0{p{y)). The notations 0{-) and Q,{-) are described 
in Section HIl 

Proof of Theorem [7] [Upper bounds;] We start by establishing 
the upper bounds on the effective resistance for graphs that are 
dense in R'^. Throughout the proof of the upper bounds, we 
will use R'^y{G), for any graph Q, to denote the effective 
resistance between nodes u and v in the electrical network 
{Q, -Pmax) with every edge of Q having a generalized resistance 
of Pmax- From the Electrical Analogy theorem and Monotonic- 
ity Law (Theorems m and ID, we get 

E«,o < KliG). 

To establish an upper bound on I]„ o, we will now establish an 
upper bound on the resistance R'^oiG)- To this effect, suppose 
that / is a dense drawing of G in M''. From dense embedding 
Theorem |3] we conclude that there exists a positive integer h 
such that the d-D lattice Z^ can be embedded in the /i-fuzz 
of Q. Moreover, Theorem [3] tells us that there exists Uz,Oz € 
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T^Zd, a positive constant c, and an embedding i] : l/z^ ^ 
of Zd into C/C*), such that 

dg{u,r]{uz)) < dg{o,r]{o^)) < c (11) 

dza{uz,0:,) < 4:d+ —df{u,o), (12) 
7 

where 7 is the maximum uncovered diameter of the /-drawing 
of Q. Note that ri{uz),ri{oz) G 1^ ■ Consider the electrical 
network {G^'^\ Pmavi) formed by assigning to every edge of 
G'^'*^ a resistance of Pmax- From the triangle inequality for 
effective resistances (Lemma O, 

+ ^;(o.),o(e''^)- (13) 

For any two nodes u,v e 1/", application of the triangle 
inequality Lemma [3] to successive nodes on the shortest path 
joining u and v gives us R°u,v{Q^^^) < '^gC') w)^'max < 
f )Pniax- Using this bound in ( fTST l, and by using ( fTTI ). 
we conclude that 

< l^'max + (14) 

Since CJ^''' D Z^, from Rayleigh's Monotonicity Law (Theo- 
rem |5]l, we obtain 

When Q is dense in, say, in K^, we have from Lemma |5] that 

i?°*f,<,jZ2) = e(iogdz,K,o.)), 

which implies 

^;(n.),,(0.)(e*''^)=0(l0grfz.(".,0.)). 

Combining this with (fTZt and ( fT4] i. we get 

i?f„(C/(''))=0(logd/(u,o)). 

From Lemma |4] we know that the effective resistance in Q and 
its /i-fuzz is of the same order, so that 

from which the desired result follows: 

S«,o<<o(e)-O(logd/Ko)). 

The statements of the upper bounds for 1 and 3-dimensions 
can be proved similarly. This concludes the proof of the upper 
bounds in Theorem [T] 

[Lower bounds:] Now we establish the lower bounds on the 
BLUE error covariance S„ o in a sparse graph. Throughout the 
proof of the lower bounds, for a graph Q, we will use Rlfi,{G) 
to denote the effective resistance between nodes u and v in 
the electrical network (CJ, Pmin) with every edge of Q having 
a generalized resistance of Pmin- From the Electrical Analogy 
and Rayleigh's Monotonicity Law (Theorems |4] and |5]l, we get 

S«,o > RtoiQ)- (15) 

Therefore, to establish a lower bound on we proceed 

by establishing a lower bound on the resistance R'^oi^)- 



Since Q is sparse in W^, it follows from Theorem |2] that 
there exists a positive integer h, such that Q C Z^^''''. Let 
r] : 1/' —f l/za be the embedding of Q into 7i'j^\ Consider 
the generalized electrical network (Zj^'*'' , Pmin) formed by 
assigning a generalized resistance of Pmin to every edge of 
zIj''"*. From Rayleigh's monotonicity law, we get 

KliQ)>KloA^^d^), (16) 

where Uz — r]{u),Oz — rj{o) refer to the nodes in Z^'*-' that 
correspond to the nodes u, o in Q. When the graph is sparse 
in, say, M^, it follows from (fTSI l and Lemma |5] that 

Rto{Q)^^{^Ogdz,{Uz.Oz)) 

= n{\ogdf{u,o)) , 

where the second statement follows from ^ in Theo- 
rem 12] Combining the above with (fTsT i, we get Eu.o = 
57(log o)), which proves the lower bound for graph that 
are sparse in R^. The statements for the lower bounds graphs 
that are sparse in or M"^ can be proved in an analogous 
manner This concludes the proof of the theorem. 

VI. Checking Denseness and Sparseness 

To show that a graph is dense (or sparse) in a particular 
dimension, one has to find a drawing in that dimension with 
the appropriate properties. For sensor networks, sometimes the 
natural drawing of a deployed network is sufficient for this 
purpose. By the natural drawing of a sensor network we mean 
the mapping from the nodes to their physical locations in the 
Euclidean space in which they are deployed. We can use this 
natural drawing to construct the following examples of dense 
and sparse graphs. 

Proposition 1: 1) Deploy a countable number of nodes 
in so that the maximum uncovered diameter 7 of its 
natural drawing is finite, and allow every pair of nodes 
whose Euclidean distance is no larger than 27 to have 
an edge between them. The resulting graph is weakly 
connected and dense in W''. Such a graph is also sparse 
in if the nodes are placed such that every finite volume 
in contains a finite number of nodes. 
2) Consider an initial deployment of nodes on a square 
lattice in R^, for which a fraction of the nodes has 
subsequently failed. Suppose that the number of nodes 
that failed in any given region is bounded by a linear 
function of the area of the region, i.e., that there exist 
constants a and (5 such that, for every region of area A 
the number of nodes that failed in that region is no larger 
than a A + /3. Assuming that a < ^p^^j , there will be 
an infinite connected component among the remaining 
nodes, which is dense and sparse in 2-D. □ 

The proof of the proposition above is provided in Ap- 
pendix u 

The first example in the proposition is that of a geometric 
graph that is obtained by placing a number of nodes in a region 
and specifying a range such that a pair of nodes have an edge 
between them if and only if the Euclidean distance between 
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them is no more than the given range. The second example 
refers to a network in which some of the initially deployed 
nodes have failed, with the stipulation that in large areas, no 
more than a certain fraction of the node may fail. For example, 
/? = 5 and a = 0.04 satisfies the stated conditions. It can be 
shown that [3 ~ 5 and a — 0.04 means that in areas larger 
than 10 X 10, at most 4% of the nodes may fail. 

To show that a graph is not dense (or not sparse) in a 
particular dimension is harder since one has to show that 
no drawing with the required properties exists. Typically, this 
can be done by showing that the existence of a dense (or 
sparse) drawing leads to a contradiction. An application of 
this technique leads to the following result. 

Lemma 6: 1) The c?-dimensional lattice is not sparse 
in K- for every d < d, and it is not dense in M.'^ for every 
d> d. 

2) A regular-degre^ infinite tree is not dense or sparse in 
any dimension. □ 

The first statement of the lemma is provided in Appendix J] 
The proof of the second statement is not provided since the 
method of the proof is similar. 

We are now ready to prove Lemma |2l 

Proof of Lemma |2] To prove the result by contradiction, sup- 
pose that a graph Q is dense in as well as sparse in M"* , 
where d' < d. It follows from Theorems [3] and |2] that there 
exist positive integers £,p such that C G^^'' and Q C zj^^?'. 
It is straightforward to verify the following facts: 

1) for every pair of graphs G, Q that do not have any parallel 
edges, Q (z Q ^ t/*^'-* C Q^^^ for every positive integer I. 

2) for an arbitrary graph Q without parallel edges, and two 
positive integers i,p, we have (C/(p))(^) — gip^) _ 

It follows that Zrf C Z^^F\ which means, from sparse 
embedding Theorem |2] that a d-dimensional lattice is sparse 
in K*^ . This is a contradiction because of Lemma |6] which 
completes the proof. 

VII. Summary and Future Work 

In a large number of sensor and ad-hoc network applica- 
tions, a number of node variables need to be estimated from 
measurements of the noisy differences between them. This 
estimation problem is naturally posed in terms of a graph. 

We established a classification of graphs, namely, dense or 
sparse in M'^, 1 < d < 3, that determines how the optimal lin- 
ear unbiased estimation error of a node grows with its distance 
from the reference node. The notion of denseness/sparseness 
introduced in this paper is distinct from the usual notion based 
on the average degree. In fact, we illustrated through examples 
that node degree is a poor measure of how the estimation error 
scales with distance. 

The bounds and the associated graph classification derived 
here can be used in performance analysis, design and deploy- 
ment of large networks. For example, if a sensor network 

^A graph is called regular-degree if the degree of every node in the graph 
is the same. 



is sparse in M, then we know that the estimation error of a 
node will grow linearly with its distance from a reference. A 
large number of reference nodes will thus be needed for large 
networks that are sparse in M. On the other hand, if one has 
control over the network deployment, then one should strive 
to obtain a network that is dense in with d as large as 
possible. In the ideal case of c? = 3, with a single reference 
node one can get bounded estimation error regardless of how 
large the network is. 

There are several avenues for future research. The scaling 
laws described in this paper were derived for infinite mea- 
surement graphs. This is justified by the fact that the BLUE 
covariance Y^^.o of a node u in an infinite graph is very close 
to the obtained in a large finite subgraph that contains the 
nodes u and o sufficiently inside it |7l]. However, to gain 
a better understanding of the "boundary" effects that can 
occur in finite graphs, an interesting research direction is to 
determine how large the BLUE error covariance can be as 
a function of the size of the graph, for nodes that are close 
to the edge of the graph. A connection between the notions 
introduced in this paper and those in coarse geometry might 
be useful in this regard. It can be shown that a graph Q 
that is both sparse and dense in is coarsely equivalent 
to W^, which intuitively means that Q and are the same 
in their large scale structure (see ll23ll for a precise definition 
of coarse equivalence). Certain coarse geometric notions that 
were originally defined for infinite graphs have been extended 
to finite graphs (see [24]). This connection between coarse 
geometry and denseness/sparseness might provide a way to 
extend the techniques used in this paper to finite graphs. 

Although the dense and sparse classification does allow 
randomness in the structure of the graph, the effect of such 
randomness on the scaling laws for the error is not explicitly 
accounted for in the present work. A useful research direction 
would be to characterize the estimation error covariances 
in graphs with random structure, such as random geometric 
graphs [25]. Another interesting avenue for future research 
is the investigation of estimation error growth in scale-free 
networks that do not satisfy the bounded degree assumption. 
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Appendix I 
Technical Proofs 

Proof of Lemma \T\ We prove that [T] implies |2] by contradic- 
tion. Assuming that |2] does not hold, we have that 

Va > V/3 > 3u,v e 'U such that dg{u,v) > adf{u,v) 

or equivalently 

Va > V/3 > 3u,v£ 1^ 

such that 

df{u,v) _ (3 
dg{u,v) a adg{u,v) 
This means that for a given a > 0, /3 > 0, the set 

dfiu,v) 

■ dg{u, v) 

contains at least the element 



(44^^ : u,w e 0^ and dg{u,v) > /?) 
lda{u,v) J 



df{u,v) _ 



1 

< - 



dg{u,v) a adg{u,v) a 



and therefore 

inf (44-^ ■■u.v and dg{u,v) > (iX <-. 
ldg(u,v) J a 

Making P ^ oo we obtain that 

p = lim inf (44~^ : u,v e and dg{u, v)> (i\ <—. 
p^oo ldg{u,v) J a 

But since a can be arbitrarily large, the above actually implies 
that (0 = 0, which contradicts [T] 

To prove that |2] implies [T] we note that when |2] holds, we 
conclude that for every pair of nodes u,v £ 'V , for which 

dg{u,v) > n, we have that 



df{u,v) ^ 1_ 



fi 



1 13 

> — -. Vu^v e I'. 



dg{u,v) ~ a dg{u,v) ~ a n' 
Therefore, 

. , r df(u, v) , . , / N 1 1 /? 

iiif i -ri \ : u,v £ 1^ and dg{u,v) > n} > . 

ldg[u,v) Jan 

As n oo, the left-hand side converges to p and the right- 
hand side converges to ^ > 0, from which [T] follows. 

Proof of Theorem \3] In this proof, we will denote by g : 
I'za K'^ the natural drawing of the lattice Z^. 
(^) We have to prove that if Q is dense in M.'^, conditions 
(i) and (ii) are satisfied. Since Q is dense in M'', there is a 
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drawing function f : 1/ ^ M."^ so that the /-drawing of Q has 
a 7 < cxD and p > 0. Define a new drawing /' : 1^ — > R'' as 

7 

so that the maximum uncovered diameter 7' of the /' drawing 
of G is 1. Note that /' is still a dense drawing of Q. Now we 
superimpose the natural g-drawing of Zd on the /'-drawing of 
Q, and draw open balls of diameter 1 centered at the natural 
drawing g{uz) of every lattice node, denoted by B{g{uz), 
Figure |5] shows an example in R^. Since 7' = 1, it follows 
from the definition of denseness that in every one of those 
balls, there is at least one node u E 1^. To construct the 




Fig. 5. Superimposing a 2-dimensional lattice (gray/brown) on a 2- 
dimensional dense grapli (black). 



embedding, we associate each node of the lattice to a node of 
Q whose drawing appears inside the ball centered around the 
lattice node. This defines an injective function 7/ : l^z^ 1^. 
Consider two nodes of the lattice Uz,Vz € l^z^ that have 
an edge between them. Let u := ri{uz),v := rj{vz)- Since 
f'{u) and f'{v) belong to adjacent balls of unit diameter (see 
Figure |5]l, 

df,{u,v) = \\r{u)- f'{v)\\<2. 

Since /' is a dense drawing in R'' with 7' = 1, it follows from 
Lemma [T] that dg{u^ u) < 2a + (3, for some positive constants 
a and /3. Define h := \2a + (3~\. Then u and v will have an 
edge between them in the /i-fuzz Q'^'^K So ^'^''^ D Z^, and we 
have the desired result that denseness implies (i). 

To show that denseness implies (ii), first note that if u G 
vi'^Za)^ then (ii) is trivially true (choose u :— u), so only 
nodes in 1/" \ rj{l/xi) are interesting. For every u E 'V , find 
Uz E I^Zd as the node in the lattice such that the ball of 
unit diameter drawn around Uz is closest to u. That is, find 
Uz € I^Zd such that 

Uz = arg min dist (/'(w), S(.g(wl), 1/2)) (17) 

where dist (a;, A) between a point a; G M"* and a set A C 
is defined as 

distfa;, A) = inf lla; — y\\. 
yeA 



There are only 2'^ balls one needs to check to determine 
the minimum in ([IT), so Uz exists, though it may not be 
unique. If there are multiple minima in ([IT), pick any one. 
This procedure defines an onto map £^ : 'U "^Za- Let 
1] : I^Zd ^ 1^ be the embedding of into Q''^'> as described 
earlier in this proof. Define ifj : 'V' ^ 1^ as ip :— {r] o . We 
will now show that, for every u E 'P' , the node ijj{u) E I' , 
which has a corresponding node in the lattice, is within a 
uniformly bounded graphical distance of u. Since f'{u) either 
lies in the ball centered at g{uz) or in the gaps between 
that ball and the neighboring balls, \\f'{u) — g{uz)\\ < \fd- 
Therefore, 



ds,{u,i;{u)) < Wfiu) - g{uz)\\ + \\g(,Uz) - /'(VX"))II 



(18) 



where we have used the fact that f'{'ijj{u)) E B{g{uz), 
From Lemma [T| and the denseness of the /-drawing of Q, we 



get 



dg{u, ip{u)) < adf{u, ii{u)) + (3 
— a^fdfi (u, + (3 

3 



Define 



< -a7\/d + /3. 



c:= [-a7yd + /3], 



(19) 



which is a constant independent of u and v. Then for every 
u E 'V , there exists a. u := ip{u) E 'ri{l/zd) ^ such that 
dg{u, u) < c, which is the desired condition (ii). 

(<J=) We have to prove that if (i) and (ii) are satisfied, then Q is 
dense in R''. We will construct a drawing f of Q in M.'^ with the 
following procedure and then prove that it is a dense drawing. 
Since Z C G^'^\ there is an injective map 77 : I^Zd ~^ 
such that ri{'V'za) C 1^. Pick a node u in 'P' that has not been 
drawn yet. By (ii), there exists a positive constant c and a 
node Uz G 'P'za such that u :— ri{uz) E V and dg(u^ u) < c. 
If u has not been drawn yet, then draw it the location of its 
corresponding lattice node, i.e.. 



f{u) = g{uz). 



(20) 



A little thought will reveal that if u has been drawn already, 
as long as the drawing procedure outlined so far is followed, 
it must have been drawn on the lattice location g{uz), so (|20l ) 
holds. Once u is drawn, we draw u in the following way. In 
case u — u, drawing of u is determined by the drawing of 
u. If u ^ u, draw u by choosing a random location inside 
an open ball of diameter 1 with the center at f{u). To show 
that a drawing obtained this way is dense, first note that the 
largest uncovered diameter 7 < 2 since a subset of the nodes 
of occupy the lattice node positions. Pick any two nodes 
u,v E V . Again, from (ii), we know that there exists u, w G 
vi'^Za) C 1^ such that dg{u,u) < c and dg{v,v) < c for 
some positive constant c. Therefore 

dg{u, v) < dg{u, u) + dg{u, v) + dg{v, v) 
<2c + hdg(h){u,v) 
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Since C G^^K 

= hiuz) - g{vz)\\i 
where || • |li denotes the vector 1-norm, 

< Vd\\giu^) - g{vz)\\ 

= Vd\\f{u) - f{v))\\ (from (|20|) 

= Vddfiu, v). 

Because of the way the drawing / is constructed, we have 
df{u, u) < 1, which implies df{u, u) < df{u, u) + df{u, v) + 
df{v, v) — df{u, v) + 2. So we have 

dg{u,v) < 2c+ hVd{df{u,v) + 2) 

= 2(c + hVd) + hVd df{u, v). 

From Lemma [T] we see that the asymptotic distance ratio p > 
for the /-drawing of Q, which establishes that / is a dense 
drawing of Q in M'*. It follows that Q is dense in M''. 

To prove the relationship ^ for any dense drawing /, consider 
again the scaled drawing /' defined as /' — f /j, so that the 
maximum uncovered diameter of /' is 1. Since Q is dense 
in R'', Zd can be embedded in Q^'^'^ with an embedding 77 : 
I^Zd ■ We choose the embedding rj as described in the 
first part of the proof. For every u 'P' , call Uz ■= £,{u), 
where ^ : 1^ ^ 'V'za was defined earlier in this proof for 
the /' dense drawing of Q. Now consider two arbitrary nodes 
u,v G 'V' and let Uz £,{u), Vz ■= £,{v) (see Figure |6]l. It 
was shown earlier in this proof that for every pair of nodes 
u,v £ 1/ , we have dg{u,r]{uz)) < c and dg{v,r]{vz)) < c, 
where c is defined in ( fT9] l. 




7,^ = 1 



1' Vz 



V 



Fig. 6. Natural drawing of tlie 2-D lattice (gray/brown) superimposed on 
the /' drawing of Q. Edges are not shown to prevent clutter. In this example, 
u = u but V ^ V. 



Now, 



and 



dzaiuz,Vz) = \\giuz) - 9{Vz)\\l 

< Vd\\giuz) - g{vz)\\, 



Muz) g{vz)\\ < \\9{uz) - + - f'{u)\\ + 

\\f'{u)^f'{v)\\ + \\f'{v)-f'{v)\\ + 
\\nv)-g{vz)\\. 



Fig. 7. 2D lattice with certain fraction of nodes removed to simulate node 
failure. Nodes are shown as dots, except nodes that belong to disconnected 
components created due to failed nodes, which are shown as circles. Edges 
are not shown to avoid clutter The superimposed grid has cell size equal to 
L X L, cell boundaries being shown in solid lines. 



We know that \\g{uz) - /'(w)|| < 5 < ^ since f'{u) e 
B{g{uz), i), and \\f'{u) - f {u)\\ < fVrf from dH. Using 
these in the above, we get 

\\g{uz) - g{vz)\\ < 4Vd + df,{u,v), 



^dz^iuz,Vz) < M 
which is the desired result. 



Vd 



df{u,v) 



Proof of Proposition\T\ Proving that the graph in the first 
example has the required properties is straightforward and is 
therefore omitted. For the second example created from a 2-D 
lattice after removal of a certain fraction of its nodes, we first 
have to prove there is a giant connected component. First of 
aU, a < implies there exists a positive integer L such 

that aL'^+P < L — 1. Pick such a L and superimpose a square 
grid on the natural drawing of the lattice Z2, with each side 
of the grid having length L; as shown in Figure [T] Pick an 
arbitrary cell C. The nodes in C that have failed may divide 
the remaining nodes inside the cell into multiple connected 
components. Consider the largest connected subgraph formed 
by the remaining nodes inside the cell, and call it Qc = 
{'Uc, "Ec)- The edges in Ec are incident on only those nodes 
that are entirely in T^c- There is at least one node on each of 
the four boundaries of C that is part of Qc- If not, all the L 
nodes on at least one side must either fail or be removed from 
the largest component due to separation by the failed nodes. 
In both case, the number of nodes that must fail in a region of 
area is at least L, which means aL^ + (3 > L. However, 
this violates the condition that aL^ + P < L — 1. Moreover, 
each of these boundary nodes are connected to one another by 
a path that lies entirely in Qc. For, if not, there must be a "fault 
line" created by failed nodes that divides the cell into a left and 
a right (or, top and bottom) half, which requires the removal 
of L nodes. By the argument present above, such a situation 
cannot occur. Therefore, the largest connected component in 
each cell contains at least one node from each of its four sides. 
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Now consider two adjacent cells Ci and C2, sharing a common 
side of length L. Pick a node 6 that belongs to Qci and 
lies on the common boundary between Ci and C2, which 
exists by the arguments above. This node h also belongs 
to Gc-i^ since otherwise this node has to be separated from 
the nodes in by a "fault line" of failed/removed nodes, 
which was impossible. Therefore, the graphs Qc^ and 
are connected. Moreover, a conservative upper bound on the 
graphical distance between two arbitrary nodes in Qc^ and 
Qc2 is 8L. Continuing this argument, we see that the largest 
connected component in each cell is connected to those in 
adjacent cells, which proves that there is a giant connected 
component formed by the nodes after removal of the failed 
nodes from Z2. We call this giant component Q. To prove that 
Q is dense in 2-D, first consider the natural drawing of Q that 
is obtained from that of Z2 after removal of the failed nodes. 
The largest uncovered diameter 7 in this natural drawing of Q 
is trivially finite. To prove that p > 0, pick two arbitrary nodes 
u and V in Q. The minimum number of adjacent cells of size 
L X L one has to go through to go from u to w is no greater 
than [ ^'^^^"'"'' ], where df{-,-) is the Euclidean distance in 
the natural drawing of the lattice Z2. Therefore, the graphical 
distance between u and v satisfies 

dg{u,v) < 8L[^^^^4^1 < 20df{u,v). 

Lemma [T] now implies that p > 0, which proves denseness 
of Q. Sparseness of Q follows trivially from the sparseness of 
Z2. 



Proof of Lemma |6] We only provide the proof that the 2- 
dimensional lattice is not sparse in M and is not dense in M?. 
The general case for arbitrary dimensions is analogous. 
To prove by contradiction the lack of denseness, assume that 
there exists a dense drawing / of Z2 in M"^, with associated 
7 < cx) and p > {). Fix the origin of at f(u) for an 
arbitrary node u in the lattice Z2. For an arbitrary D > Q, the 
volume of the sphere in 'E? centered at the origin with diameter 
D, denoted by ^3^(0, £>), is rt{D'^). Therefore the number of 
nodes of Z2 drawn inside B^[0,D) is Vl{{^f) = n{D^). It 
is straightforward to show that for any set of n distinct nodes 
in the lattice Z2, the maximum graphical distance between 
any two nodes in the set is f2(-yn). Therefore the maximum 
graphical distance between the nodes in S'^(0, D) is f2(Z?2 ). 
The maximum Euclidean distance between any two nodes 
drawn inside the sphere S'^(0, D) under the /-drawing is at 
most D, and since / is a dense drawing, it follows from 
Lemma [T] that for every pair of nodes u, v in Z2 such 
that f{u)J{v) e B^iO.D), we have dG{u,v) < aD + b. 
Therefore, the maximum graphical distance between pairs of 
nodes whose drawing falls inside B^{0,D) is 0{D), as well 
as il{D^), which is a contradiction for sufficiently large D. 
Hence no dense drawing of Z2 in R'^ is possible. 

To show Z2 is not sparse in R, assume that there exists a 
civilized drawing of Z2 in R with s > and r < 00, where 
r and s are constants. Consider a subgraph Z2(„) of Z2 that 
consists of all nodes within a Euclidean distance n from the 



origin. The total number of nodes in this finite subgraph is 
r2(n^). The length of the interval, L, in which the nodes of this 
subgraph are located in the sparse 1-d drawing of Z2 is clearly 
L = f7(sn^). Since the maximum graphical distance between 
any two nodes in the subgraph Z2(„) is n by construction, the 
maximum connected range in the 1-d drawing must be at least 
r > ^ = il{sn). Since this must be true for every n, r cannot 
be a finite constant. Thus, no civilized drawing of Z2 in R 
exists. 

Appendix II 

Relative position measurements from noisy range 

AND angle measurements 

Here we show how to convert noisy measurements of range 
and angle between a pair of nodes into a measurement of 
relative position between them such that the error in the 
measurement is zero mean. 

Let the position of two nodes u and v lying on a plane (in a 
in a 2D Cartesian reference frame) be denoted by and Xy, 
respectively. Let the relative position between them (i.e., x„ — 
Xy) be [Ax, Ay]^, and the distance and angle between them 
(in the same reference frame) be r and 9. This implies Aa: = 
r cos 6 and Ay = r sin 9. Let r and 9 be noisy measurements 
of the distance and angle: r = r + Sr, 6 = 6 + 69, where 6r 
and 69 are random errors in range and angle measurements. 



The following is an unbiased measurement of Xy 



Xfj . 



Cu,v 



Ax 


1 


rcosO 


Ay_ 


c 


rsinfl 



where c := E[cos((5e')], (21) 



under the assumption that (i) 5r and 56 are independent, (ii) 
5r is zero mean, and (iii) E[sin(5^?] — 0. 

Before we show why the measurement error (u,v — {xu — 
Xy) is zero mean, we note that the measurement C„ „ can be 
obtained in practice from range and angle measurements and 
knowledge of c. The value of c can be obtained (estimated) 
a-priori while characterizing the AoA sensor used to gather 
angle measurements. For example, if 69 ~ N{Q, cr^), then 
c ~ g-o.5a- (yyijiQij Q^jj verified by symbolic integration in 
MATLAB). 

To show that the measurement error (u,v ~ [Ax, Ay]-^ is 
indeed zero mean, we expand the expression for Ax to obtain 

1 



Ax 



(Axcos((56') - Aysm{5e) + 5rcos{9 + 69)) (22) 



We decompose cos 69 as 



■ C50, 



(23) 



where E[c56/] = by construction. From (|22] |. we get 

Ax = Ax + - {cseAx — Ay sin{69) + 6r cos( 



69)), 
(24) 



and similarly. 

Ay = Ay- 



- (c^e Ay + Ax sin < 
c 



6rsm{9 + 69)), (25) 
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which yields 



Ax 




'Ax 










Ay 


+ 





where 

fix ■= \ (cseAx - Ay sin SO + 6r cos(9 + 69)) (26) 
c 

ey ■■= -r {cseAy + Ax sin 60 + 6r sm{e + 66)) . (27) 

The random variables and ey are zero mean due to the 
assumptions on 69 and 6r, which verifies the claim that the 
measurement error Cu,v — {xu — x^) is indeed zero mean. 
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